Journal of Nonlinear Mathematical Physics Volume *, Number * f 200*1, [T hill 



Article 



Lower bounds for the spinless Salpeter equation 

Fabian Brau 1 

Service de Physique Generate et de Physique des Particules Elementaires, Groupe de 
Physique Nucleaire Theorique, Universite de Mons-Hainaut, B-7000 Mons, Belgique 
E-mail: fabian.brau@umh.ac.be 

-rf ', Received Month *, 200*; Revised Month *, 200*; Accepted Month *, 200* 

o i 

Abstract 

<N 

We obtain lower bounds on the ground state energy, in one and three dimensions, for 
the spinless Salpeter equation (Schrodinger equation with a relativistic kinetic energy 
operator) applicable to potentials for which the attractive parts are in L p (W l ) for some 
p > n (n = 1 or 3). An extension to confining potentials, which are not in L p (M. n ), is 
also presented. 



(N 



> 

G\ ■ 1 Introduction 

O ' 
O 



The spinless Salpeter equation is a simple relativistic version of the Schrodinger equation 
which can be obtained, with some approximations, from the Bethe-Salpeter equation ^H]. 
A covariant description of bound states of two particles is achieved with the Bethe-Salpeter 
equation pQ. This equation reduces to the spinless Salpeter equation [2] when the following 
i-£h ■ approximations are performed: 

• elimination of any dependences on timelike variables (which leads to the Salpeter 
equation [H]). 

• any references to the spin degrees of freedom of particles are neglected as well as 
negative energy solutions. 



X 



The spinless Salpeter equation takes the form (H = c = 1) 



ay/p 2 + m 2 + V(r) *(r) = M *(r), (1.1) 

where m is the mass of the particle and M is the mass of the eigenstate (M = am + E, E 
is the binding energy). The parameter a equals 1 for a one-particle problem and equals 
2 for a two-particle (identical) problem (in the center of mass frame). We restrict our 
attention to interactions which are introduced in the free equation through the substitution 
M — > M — V(r), where V(r) is the time component of a relativistic four- vector. The 
interaction could also, in principle, be introduced through the substitution p — > p — A(r), 
where A(r) is the spatial component of a relativistic four-vector. However we do not 
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consider this kind of potentials since the derivation of the spinless Salpeter equation from 
the Bethe-Salpeter equation leads to A(r) = 0. Equation (|1.1|) is generally used when 
kinetic relativistic effects cannot be neglected and when the particles under consideration 
are bosons or when the spin of the particles is neglected or is only taken into account 
via spin-dependent interactions. Despite its apparent complexity, this equation is often 
preferred to the Klein-Gordon equation. The equation (|l.lj) appears, for example, in 
mesons and baryons spectroscopy in the context of potential models 1 (see for example 

H El)- 

Conversely to the Schrddinger equation, for which a fairly large number of results, 
giving limits on the number of bound states or limits on the values of energy levels, can 
be found in the literature (see for example [111111231111111123121113111111311111131301 
I3T1 132*1 133} E3 CHH E3 ) j there are few general rigorous results concerning the spinless 
Salpeter equation due to the pseudo-differential nature of the kinetic energy operator. 
Most of these results have been obtained for a Coulomb potential (for example, upper 
and lower bounds on energy levels) H EH El El- Recently upper and lower limits 
on energy levels have been obtained for some other particular interactions |14l 1151 fTHl I17j . 
Some more general results also exist. An upper bound on the total number of bound 
states in three dimensions and an upper limit on the number of £-wave bound states have 
been obtained in Refs. jSEHHH]- An upper limit on the critical value, gi\ of the coupling 
constant (strength), g, of the central potential, V(r) = gv(r), for which a first l-wave 
bound state appears has been derived in Ref. 10 . 

In this article we obtain lower bounds on the ground state energy, in one and three 
dimensions, for the spinless Salpeter equation applicable to potentials for which the at- 
tractive parts are in L p (M. n ) for some p > n (n = 1 or 3). These lower bounds yield lower 
limits on a second kind of critical value, g c , of the coupling constant, g, of the potential, 
V(r) = gv(r), for which the mass of the eigenstate M is vanishing, entailing that for any 
value of g greater than g c the system has an unphysical negative mass. 



2 General formulas 

In this Section, we derive some general formulas that we apply in Sections H and H to 
obtain the lower bounds in one and three dimensions. 

We begin by recalling the expression of the sharp Young's inequality obtained in 1975 
by Beckner and independently by Brascamp and Lieb |42l I43j 



d n X d n yf(x)g(x-y)h(y) 



<[C p C q C r ) n \\f\\ p \\g\\ q 



(2.1) 



where p 1 + q 1 + r 1 = 2 with p, q, r > 1 and 



/ 



<Tx|/(x)P> 



i/p 



(2.2) 



1 For a review of several aspects of the "semirelativistic" description of bound states with the spinless 
Salpeter equation see: W. Lucha and F. F. Schoberl, Semirelativistic treatment of bound states, Int. J. 
Mod. Phys. A 14 (1999), 2309-2334 and references therein. 
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The constants appearing in the inequality (|2.1j) are best possible and take the form 

(C P ) 2 =p 1/p p'~ 1/p \ (2.3) 

where p'~ l = 1 — p^ 1 . 

Now we write the spinless Salpeter equation as an integral equation 

*(x) = J <fyGW(A) [M - V(y)] *(y), (2.4) 

where A = x — y, A = |A| and G^ n \A) is the Green's function of the kinetic energy 
operator, a [p 2 + m 2 ] 1 ^ 2 , in n dimensions. For n = 3 the expression of the Green's function 
is known 41 j (see also Section |2J). For n = 1 the Green's function is obtained in Section [SJ 
The Young's inequality (|2,1[) can be used together with the expression (|2,4[) of the 
spinless Salpeter equation to obtain 

mi<\M\[CpC § Cr] n 1 1*| | f | |G (n) I |f 1 1*1 If + [C p C q C r ] n ||y-*|| p ||G( n )|U|*|| r , (2.5) 

where V _ (x) = max(0, — V(x)) is the attractive (negative) part of the potential. In order 
to obtain a formula without any reference to the wave function, we choose p = f = 2, 
q = 1 and r = 2. With the help of Holder's inequality applied to ||y _ *|| p , it is possible 
to write 

||*||1 < |M|||G (n) ||i||*j|^ + [QC^] n ||F-||^||G ( "' ) |L ||*||1 (2.6) 

Zq-2 9-1 

with q > 1. We obviously suppose for now that ||G (n) || g exists for some 1 < q < q\ . 
This question is discussed in Sections 01 and |SJ The quantity ||V r_ ||_9_ exists provided the 

9-1 

attractive part of the potential is in ^(W 1 ) for some p > q^ /(Qo — !)• 
From ()2.6|) we finally obtain the following lower limits 

|M| > HG^II^ 1 [l-Q||y-||_9_ ||G( n )||g] , (2.7a) 

where the constant takes the form 

<? = l^]' = ^(^) (5^) (^) • ("b) 

Physically M is required to be nonnegative and the quantity \M\ could be replaced by the 
mass itself, M, provided that we restrict the right-hand side of (|2.7a|) to be nonnegative. 

Equation (|2.7|) is the main result of this article. In the next Sections we apply this 
relation to obtain lower bounds in one and three dimensions. 



3 Lower bound in three dimensions 

In order to apply the relation (|2.7|) in the specific case of a three-dimensional space we 
need the relevant Green's function. The computation has already been done in The 
Green's function is obtained by computing the following integral 

G (3) (A) = 7^/^ eXP( ": P - A) . (3-D 



(2vr) 3 J ayV + m 2 
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We find the following expression 



G (3) (A) 



in 



2avr 2 A 



^i(mA), 



(3.2) 



where K v {x) is a modified Bessel's function (see for example [441 p. 374]). With the 
expression (|3.2JI for the Green's function we can compute the needed norms. We find that 



\\G^\\i = — 



am 



(3.3) 



and 



\\G®\ 



m 



2-3/5 



(47T) 



2avr 2 



dxx 2 ~ q K^xf 



m 



2-3/ q ^£ 



a 



(3.4) 



These relations, (|3.3|) and (|3.4|) . imply that we cannot consider the ultrarelativistic regime 
(m = 0) in order to have finite values for the various norms of the Green's function. 
Moreover, since the Bessel's function K\(x) behaves like x~ l for x ~ (lim^^o xK\ (x) = 
1), we have the following restriction on q in order to have a finite value of the constant C q 



P) 



1 < q < q y o' = -■ (3.5) 
The lower bound on the mass of the ground state is (we suppose M to be nonnegative) 



M > am -C^C g m 3 - 3/q \\V || (3.6) 
This lower limit is nontrivial only if ||^ _ ||_2_ exists. The restriction (|3.5|) implies that 

9-1 

the negative part of the potential must be in L P (M 3 ) for some p > 3. Note that the limit 
q I 1 of the inequality (|3.6|) gives M > am — ||y~||oo which is trivially true. 

The inequality (|3.6|) obviously yields the following lower bound on the ground state 
binding energy 

E > -C^Cqm^^WV'W^. (3.7) 

This lower limit implies that the binding energy does not increase, in modulus, faster than 
linearly with the mass of the particle and with the strength of the potential (for the class 
of potentials we consider). 

We propose now to test the lower bound 1)3.6}) with some typical potentials. For sim- 
plicity we consider only central potentials. We use three potentials for these tests: an 
exponential potential 

V(r) = -gR^expi-r/R); (3.8) 



an exponential multiplied by a power 
V(r) = —gR~ 2 rexp(—r/R); 



(3.9) 
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and a singular potential 



V(r) 



g(rR) 



1/2 



exp(— r/R). 



(3.10) 



Note that the lower bound ()3.6|) does not apply to the Yukawa potential, for example, 
since it behaves as r _1 near the origin and thus the negative part of this potential is not 
in L p (M 3 ) for any p > 3. 

The mass of the ground state depends on two dimensionless parameters, g and (3 = mR. 
It is thus rather cumbersome to compare the lower bound (|3.6|) to the exact result. We 
propose an indirect test of this lower bound by computing a lower bound on the critical 
value, g c , of the potential, V(r) = gv(r), for which the mass M of the ground state is 
vanishing, entailing that for any value of g greater than g c the system has an unphysical 
negative mass. Prom (|3.6|) we obtain the lower limit 



This inequality depends only on one parameter, j3, and can then be tested more easily. 

In Figure Q we present a comparison between the exact values of the critical coupling 
constant, g c , obtained numerically, and the lower limit (|3. llf) as a function of /3 = mR for 
the three potentials (|3.8|) - (|3.1U|) . Note that for these tests we choose to solve the spinless 
Salpeter equation with two identical particles, a = 2, only for numerical convenience (not 
to modify numerical codes). Clearly the lower bound (|3.11|) is quite cogent and follow 
closely the exact behavior for all values of f3. The accuracy is, however, less good for the 
singular potential. 

4 Extension to confining potentials 

In this Section we consider confining potentials, such as the harmonic potential, for which 
the norm 1 1 [ | _s_ does not exist. Then the lower bound ()3.6|) does not apply without 

<3-l 

suitable modifications. This kind of potential is used, for example, in mesons and baryons 
spectroscopy, in the framework of potential models, to describe quark confinement together 
with the spinless Salpeter equation El El El • 

To obtain a lower bound for this class of potentials we consider the following truncated 
potential V(r) defined as 



g c > a C 3 g C q m 2 - 3 ^\\ 



i -l 



v 



3_ 



(3.11) 



V(r) 



V(r) for V(r) < C 
C for V{r) > C. 



(4.1) 



Since V(r) is more attractive than V(r), a lower bound on the ground state energy of 
V(r) is also a lower bound on the ground state energy of V(r). To be able to compute 
the relevant norm of the interaction we consider the potential V(r) = V(r) — C. Now the 
lower bound (|3.fi[) applies to V(r) and we have 




where My denotes the mass of the ground state of the potential V(r) and where we 
consider, for a moment, the general form of the lower bound with the modulus of My. 
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Figure 1. Comparison between the exact values of the critical coupling constant, g Cl obtained 
numerically and the lower limit l|3.11|l as a function of (3 = mR for the three potentials l|3.8() , i|3.9[) 
and l|3.1()(l . denoted respectively E, PE and S. 

The relation between the potentials, V(r) and V(r), clearly implies that 

My = M v + C. (4.3) 

However we do not have a relation between My and \M V \ and we must restrict the lower 
bound (|4.2|) to nonnegative values of My, which is secured if M > 0. This restriction 
implies a limitation on the possible values of the constant C which cannot be too large. 
We can then write 

M>M V = M V + C>M + C. (4.4) 

We have two parameters, q and C, that can be used to maximize the right-hand side 
of (|4.4|) with the restrictions 1 < q < 3/2 and M > 0. Actually the best optimization is 
obtained by searching the value q* of q which yields the greatest value C* of C such that 
M = 0. In other words we search for the value q* of q which yields the greatest value C* 
of C such that 



C\ C q m 2 - 3/q ||V"||_i_ = 1. (4.5) 
This procedure yields the lower bound, 



M > C*, 



(4.6) 
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Figure 2. Comparison between the exact values of the mass M (in GeV) of the ground state of 
the logarithmic potential i|4.7[l and the lower bound (|4.6|l as a function of j3 = mR for three values 
of g. 1 GeV has been added to the masses for g = 0.1 and 1 GeV has been subtracted from the 
masses for g — 2 (to obtain a clearer presentation). 



once C* has been found. Note that the limit q J. 1 of the inequality Q4.6JI yields M > 
am + mm(V(r)) which is trivially true. 

We propose to test this lower bound with a central logarithmic potential 

V(r) = gR- 1 \n(r/R). (4.7) 

Quigg and Rosner introduced this potential in 1977 to describe the spectrum of heavy 
mesons |45| . The values of the parameters they found was R ps 2.5 GeV" 1 and g 1.8. 
Note that these values of the parameters were obtained with the Schrodinger equation and 
not with the spinless Salpeter equation. 

In Figure |21 we present a comparison between the exact values of the mass of the ground 
state of the logarithmic potential (|4.7|) and the lower bound (|4.6|) as a function of (3 = mR 
for three values of g = 0.1, 0.5, 2. More precisely we fix R to 2.5 GeV -1 and we let m 
vary between 0.4 GeV and 4 GeV, which covers approximatively the various constituent 
(effective) masses of the quarks used in potential models (from the u and d quarks to the 
b quark). We here also choose to solve numerically the spinless Salpeter equation with 
two identical particles, a = 2. Again the accuracy of the lower bound is quite satisfactory, 
especially for small values of g. 
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5 Lower bound in one dimension 

To conclude this article we derive the lower bound applicable to one-dimensional space. 
Again, in order to apply the relation (|2.7jl . we need the relevant Green's function. It is 
obtained by computing the following integral 



G (1) (A) 



1 

2tt" 



dp 



exp(— ipA) 
a\Jp 2 + m 2 



We find the following expression for the Green's function 
1 



G (1) (A) 



an 



-K (mA), 



(5.1) 



(5.2) 



where K u (x) is again a modified Bessel's function. With the expression ()5.2j) for the 
Green's function, we obtain 



||G«||i = — 



and 



am 



m 



(5.3) 



-1/9 



2 1 /? 



air 



dxK (x) ( ' 



Hl/9 



111 



-1/9 ^1 

a 



(5.4) 



These relations, (|5.3|) and (|5.4f) . imply that we cannot consider the ultrarelativistic regime 
(to = 0) in order to have finite values for the various norms of the Green's function. The 
Bessel's function Kq(x) is characterized by a logarithmic singularity for x ~ and the 
constant C q is well defined for q > 1 (q^ = oo). 

The lower bound on the ground state energy is, where M is supposed to be nonnegative, 



M > am-C q C q m l - 1/q \\V~\\^_. 

This lower limit is nontrivial only if \\V~ 
potential must be in L P (R) for some p > 1. 



(5.5) 



| g exists. Thus the negative part of the 

q-l 
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